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We study the statistical properties of wavefunctions in a chaotic billiard that is opened up to 
the outside world. Upon increasing the openings, the billiard wavefunctions cross over from real 
to complex. Each wavefunction is characterized by a phase rigidity, which is itself a fluctuating 
quantity. We calculate the probability distribution of the phase rigidity and discuss how phase 
rigidity fluctuations cause long-range correlations of intensity and current density. We also find that 
phase rigidities for wavefunctions with different incoming wave boundary conditions are statistically 
correlated. 

PACS numbers: 05.45.Mt, 05.60.Gg, 73.23.Ad 



Microwave cavities have been used as a quantitative ex- 
perimental testing ground for theories of quantum chaos 
In quasi two-dimensional cavities, the component of 
the electric field perpendicular to the surface of the cav- 
ity satisfies a scalar Helmholtz equation that is formally 
equivalent to the Schrodinger equation. Complex field 
patterns, which model the wavefunction of an electron in 
a magnetic field, can be obtained making judicious use 
of magneto-optical effects 0, El • Alternatively, complex 
"wavefunctions" can be observed as travelling waves in 
open microwave cavities 0, Measured distributions 
of real and complex wavefunctions in microwave cavities 
with chaotic ray dynamics, where, traditionally, "com- 
plex" means that time-reversal symmetry is fully broken 
and the phase of the wavefunction has no long-range cor- 
relations, agree with a theoretical description in terms 
of a random superposition of plane waves [|| , as well as 
with random matrix theory and supersymmetric field 
theories 0- 

Recently, it has become possible to study the full real- 
to-complex crossover using microwave techniques 0, . 
The crossover regime is qualitatively different from the 
"pure" cases of real or fully complex wavefunctions. Un- 
like in the pure cases, the statistical distribution of wave- 
functions in the crossover regime depends on the choice of 
the ensemble: whether variations are taken with respect 
to the coordinate r, the frequency u>, or both. Whereas 
the theoretical work has been roughly equally divided 
between the two approaches, experiments usually need 
the additional av erag e over frequency to obtain sufficient 
statistics 0- 0- 03- El (see, however, Ref. for an excep- 
tion) . 

In general, a complex wavefunction may be written as 



V>(r) =e^(Vr(r)+^i(r)), 



(1) 



where ip T and ipi are orthogonal but need not have the 
same normalization [Tlj . The ratio of ip T and ipi is pa- 
rameterized in terms of the normalized scalar product of 
and its time-reversed, 
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The square modulus \p\ 2 is known as the "phase rigid- 
ity" of the wavefunction vb ^3] ■ Real wavefunctions have 
p = 1, whereas p = if vb is fully complex, i.e., vb T and 
tpi have the same magnitude. If the average is taken 
over the coordinate r only, whereas the frequency u> of 
the wavefunction is kept fixed, the wavefunction distri- 
bution follows by describing ipr and as random su- 
perpositions of standing waves kLhJ,hJ. The resulting 
wavefunction distribution depends parametrically on the 
phase rigidity \p\ 2 . Using a microwave billiard with a 
movable antenna, Barth and Stockmann have measured 
such a "single- wavefunction distribution" and found good 
agreement with the theory, obtaining p from an indepen- 
dent measurement [jj. It is the fact that p is different 
for each wavefunction that leads to the different results 
for averages over r only and over both r and ui. A cal- 
culation of averages with respect to frequency requires a 
theory of the probability distribution of p. Such a full 
wavefunction distribution, which needs theoretical input 
beyond the ansatz that each wavefunction -0 is a ran- 
dom superposition of plane waves, was first calculated 
by Sommers and Iida for the Pandey-Mehta Hamiltonian 
from random-matrix theory [lfij and by Fal'ko and Efe- 
tov [lfl for a disordered quantum dot in a uniform 
magnetic field. 

Fluctuations of the phase rigidity \p\ 2 have been iden- 
tified as the root cause for several striking phenomena 
in the crossover regime, such as long-range intensity cor- 
relations mf and a non-Gaussian distribution of level 
velocities 12] . Further, the existence of correlations be- 
tween phase rigidities of different wavefunctions causes 
long-range correlations between wavefunctions at differ- 
ent frequencies |18(. The experimental verification of 
these effects addresses aspects of random wavefunctions 
that have not previously been tested. The relative mag- 
nitude of the phase rigidity fluctuations is numerically 
small, leading to long-range wavefunction correlations of 
order of 10 percent or less 0, . This could explain 
why intensity distributions measured by Chung et al. 
could not distinguish between theories with and without 
phase rigidity fluctuations @. 
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In this letter, we consider wavefunctions in a billiard 
that is opened up to the outside world and calculate the 
probability distribution of phase rigidities for this case. 
Although time-reversal symmetry is not broken on the 
level of the wave equation itself, it is broken by the fact 
that one looks at a scattering state with incoming flux in 
one waveguide only Q . As we show here, random wave- 
functions in open cavities also have a fluctuating phase 
rigidity, and, hence, exhibit the same variety of phenom- 
ena as those in cavities with broken time-reversal symme- 
try, while they are much easier to generate in microwave 
experiments [5j- An additional advantage of the open- 
billiard geometry is the absence of fit parameters: The 
only parameter entering the wave- function distribution is 
the total number iV of propagating modes in the waveg- 
uides between the billiard and the outside world, which 
can be measured independently. 

Following previous works on this subject, we consider 
the parameter regime in which the frequency average is 
taken over a window Auj -C c/L -C u>, where c is the 
velocity of wave propagation and L the size of the bil- 
liard, and in which the openings occupy only a small frac- 
tion of the billiard's boundary. It is only in this regime 
that wavefunctions have a universal distribution and a 
description in terms of a random superposition of plane 
waves is appropriate. We limit ourselves to (quasi) two- 
dimensional billiards, in which the electric field perpen- 
dicular to the billiard plane is identified with the wave- 
function ip and the Poynting vector with the current den- 
sity j oc ]mip*Vip [l9|. A calculation of wavefunctions 
inside an open billiard is complementary to a transport 
study, for which one is primarily intersted in the relation 
between amplitudes of ingoing and outgoing waves in the 
waveguides attached to the billiard, not in the wavefunc- 
tion inside the cavity. Single-wavefunction statistics in 
open billiards in the universal regime was first consid- 
ered by Pnini and Shapiro Q and subsequently by Ishio 
et al. [2(1 I2H ]. Experimentally, wavefunctions in open 
billiards were investigated by Barth and Stockmann pj . 

The key to the calculation of P(p) in an open cavity 
is a relation between the scalar products of the in-cavity 
parts of scattering states •0,,, and tp„ and the Wigner- 
Smith time-delay matrix Q 22], 



(3) 



cavity 



where the scattering states have been normalized to unit 
incoming flux. Here the index p — 1,...,N labels 
the waveguide and the transverse mode from which the 
field is injected into the cavity. The time-delay matrix 
Q = —iS^dS/duj is the derivative of the scattering ma- 
trix S. In order to calculate the scalar product of 
the scattering state ipfi and its time-reversed we per- 
form a unitary transformation U that diagonalizes the 
Wigner-Smith time-delay matrix Q and rotates the scat- 



tering matrix S to the unit matrix |23| . 

S = U T U, Q = U^diag{T 1 ,...,r N )U. (4) 

The positive numbers i = 1, . . . , N, are the "proper 
delay times" , the eigenvalues of the Wigner-Smith time- 
delay matrix. Note that the incoming modes are trans- 
formed according to the unitary transformation U, while 
the outgoing modes transform according to U* , as re- 
quired by time-reversal symmetry. In the transformed 
basis, all scattering states are standing waves and, hence, 
have pjj = 1. Transforming back to the original basis, 
we find 



^.7=1 u :. 



2 T . 



E 



N 



(5) 



The joint distribution of the scattering matrix S and 
the Wigner-Smith time-delay matrix Q of a chaotic bil- 
liard is known from random-matrix theo ry l24l : The dis- 
tribution of the proper time delays n is |2j; 



JV 



3=1 



(0) 



where r av is the average delay time and 6{x) = 1 for 
x > and otherwise, whereas the unitary matrix U 
is uniformly distributed in the group of unitary N x N 
matrices. A numerical evaluation of the probability dis- 
tribution of the phase rigidity \p\ 2 is shown in Fig^for 
several values of N. For the case N = 2 of two single- 
mode waveguides, the probability distribution P(p) can 
be found in closed form, 

For a billiard coupled to the outside world via TV ^> 1 
channels P(p) becomes Gaussian, 



P(p) = £e-^' 2 / 4 . 



(8) 



This is the same functional form as the phase-rigidity dis- 
tribution for a quantum dot in a large uniform magnetic 
field E 0. 

Following Refs. 13. Efl the joint distributions of in- 
tensities and current densities away from the boundary 
of the cavity for one wavefunction ip^ can then be cal- 
culated from Berry's ansatz that ip^ can be written as a 
random superposition of plane waves 0] , 



ik-r 



(9) 
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FIG. 1: Probability distribution of the phase rigidity \p\ 
for a wavefunction in an open chaotic billiard, for different 
numbers of propagating modes connecting the billiard to the 
outside world. From bottom to top at the left end of the 
figure, curves correspond to N = 2, 3, 4, 6, 8, 10, 15, and 20. 
Inset: schematic drawing of billiard and waveguides. 



In Eq. @, all wavevectors k have the same modulus, 
while the amplitudes a M (k) are random complex num- 
bers. For a closed cavity, amplitudes of time-reversed 
plane waves are related, a M (k) = e 2ll ^a At (— k)*, where (j> 
does not depend on k. For an open cavity, no such strict 
relation exists, although some degree of correlation be- 
tween a M (k) and a /i (— k) persists in order to ensure the 



correct value of the scalar product of ip^ and 



E k Mk)l 2 ' 



Cf. Eq. 



(10) 



Taking the amplitudes corresponding to wavevectors 
pointing in different directions from identical and inde- 
pendent distributions, we see that Eq. l|lUfl implies a re- 
lation between the second moments of the amplitude dis- 
tribution, 



(a M (k)a M (-k)) = p w (|a p (k)| : 



(11) 



This, together with the normalization condition 
^ k (|a(k)| 2 ) = 1/A, where A is the area of the billiard, 
and the central limit theorem, provides sufficient infor- 
mation to determine the full distribution of the wave- 
function tp. 

As an example, we consider the joint distribution of 
the normalized intensity 7(r) = |-0(r)| 2 A and the mag- 
nitude of the normalized current density J = |j(r')|, 
j = (A/k)lmtfj*Wifj at the positions r and r' where 
k — lo/c. The single- wavefunction distribution factor- 
izes into separate probability distributions for / and J 
that each depend parametrically on the phase rigidity 



P P [I(r),J(r')} 



8.7 



(l-|p| 2 )3A 



2JV2 



(12) 



X J 



I\p\ 



1- P : 
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where Iq and Kq are Bessel functions. When both po- 
sition and frequency are varied to obtain the ensemble 
average, a further average over p is required, 

P(/(r), J(r')) = / d P P(p)P p (I(r), J(r'))- (13) 



After such average, P(I, J) no longer factorizes. The 
degree of correlation is measured through the correlator 



(/(r) 2 J(r') 2 ) c = 



-^ Var |p| 2 : 



(14) 



where (AB) C = (AB) — (A)(B) denotes the connected 
average. (Since normalization implies that (J (r)) = 1 for 
each wavefunction, correlators involving the first power of 
I factorize.) For a billiard with two single- mode waveg- 
uides, var|p| 2 = 8(148 In 2 - 128(ln2) 2 - 41)/9 « 0.078, 
cf . Eq. (0 . Similarly, we find for the correlator of inten- 
sities 



(J(r) 2 /(r') 2 ) c = var|p| 5 



(15) 



plus additional terms that describe short-range correla- 
tions. 

Thus far we have studied the distribution of a single 
scattering state in an open billiard. However, for a bil- 
liard that is coupled to the outside world via, in total, 
N propagating modes, there are N orthogonal scattering 
states at each frequency. In the remainder of this letter 
we address the question of possible correlations between 
these scattering states. 

This question can be studied using the framework of 
Ref.El which generalizes the above considerations to the 
problem of correlations between wavefunctions. Again, 
the starting point is Berry's ansatz ©, with a differ- 
ent set of amplitudes a^(k) for each scattering state 
p = l,...,N. We continue to take amplitudes a M (k) 
from identical and independent distributions for different 
directions of k, whereas we allow for correlations between 
amplitudes of time-reversed waves and between ampli- 
tudes of different scattering states. Such correlations are 
necessary, because the in-cavity parts of different scatter- 
ing states and their time-reversed states are not orthogo- 
nal, see, e.g., Eq. J2J. Hence, the second moments of the 
amplitudes a At (k) should be chosen such that 

EkM k K( k )* 



Pliv 



(E k K(k)l 2 ) 1 / 2 (E k l«,(k)l 2 ) 1 / 2 ' 

/rfr^(r)*Vy(r) 
(fdr|^(r)|2 fdr'|Vv(r')| 2 ) 1/2 ' 

Ek a M<XM- k ) 

(£j aAI (k)| 2 )i/2 (EkMk) | 2) i/2 

/rfr^ M (r)^(r) 

(/dr|^(r)| 2 /dr'|^(r')l 2 ) 1/2 ' 



(16) 



(17) 
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where, as before, the integrals are taken over the bil- 
liard only and we have chosen the normalization such 
that rinu = 1. Equations (|16fl and l|17l) then impose the 
following relations for second moments of the amplitude 
distributions: 

(a M (k)a„(k)*> = ri^Mk)! 2 ), (18) 

MkK(-k)) = /^(Mk)! 2 }. (19) 

Repeating the same arguments as those leading to Eq. 
(JSJl, we find that and p M „ can be expressed in terms 
of eigenvectors and eigenvalues of the time-delay matrix, 



1 fiiy 



P^ 



{Y,i\u j »\H j Y,i\Ui V ?Ti) 1 /> 



E 7 ' Uj/j.Uj I/ Tj 



(E.I^mIS-EJ^I^) 1 / 2 ' 



(20) 



The full distribution of the complex numbers n^ v and 
Ppu then follow from the known distributions of the 
N x N unitary matrix U and the proper time delays 
Tj , j = 1, . . . , N . A simple expression is obtained in the 
limit N 3> 1, when n^ u and p^ u acquire a Gaussian dis- 
tribution, with zero mean and with variance given by 



1 

N' 
2 



(n^n aT ) = —5^ T 5 va if P v -, 



(n^Par) = (PiiuPra) = 0. 



(21) 



Short range correlations between different scattering 
modes arise from the fact that p^ and are nonzero 
for p ^ v. These correlations exist both if statistics is 
taken as a function of position only and if the ensem- 
ble also involves a frequency average. For example, for 
the second moment of the intensity and current density 
distributions, we find from Eq. @ 

{I,{t)I„{t'))c = (|n^| 2 + |^| 2 )J (fc|r-r'|), (22) 

(jM,«( r )i^( r ')) = 4<MI?Vl 2 + l/Vl 2 ) J o(fc|r - r'|), 

with a, j3 = x,y. For the case N = 2 of a billiard with 
two single-mode waveguides one has ( | /O12 1 2 ) = (64 In 2 — 
37)/15 w 0.49 and (K 2 | 2 ) = (26 - 321n2)/15 0.25. 
Long-range correlations arise from the fluctuations of the 
"scalar products" n^ v and p^ v and exist only if the en- 
semble involves a frequency average. The lowest moment 
with long-range correlations is 

(I»(rfl v (r'f) c = -2{I,{vfj u {v'f) c (23) 

= (Ip w | 2 I/wI 2 > - (\Pw?){WA 2 )- 

For N = 2 one has (|pn| 2 |p 22 | 2 ) - (|pn| 2 ) {\p 2 2?) = 
8(5792 In 2 - 4480(ln2) 2 - 1861)/315 w 0.032. 



In conclusion, we have calculated the wavefunction dis- 
tribution for wavefunctions in an open chaotic billiard 
for the case that the ensemble average involves both an 
average over frequency and position. Fluctuations and 
correlations of the phase rigidities lead to long range cor- 
relations between intensities and current densities. Our 
results are relevant for a fit-parameter free measurement 
of the real-to-complex wavefunction crossover. 
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